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ABSTRACT 

We report on multifrequency observations performed during 2012 December-2013 
August of the first narrow-line Seyfert 1 galaxy detected in 7 rays, PMN J0948+0022 
(z = 0.5846). A 7-ray flare was observed by the Large Area Telescope on board Fermi 
during 2012 Decenrber-2013 January, reaching a daily peak flux in the 0.1-100 GeV 
energy range of (155 ± 31)xl0~ 8 ph cnr -2 s -1 on 2013 January 1, corresponding 
to an apparent isotropic luminosity of ~1.5xl0 48 erg s _1 . The 7-ray flaring period 
triggered Swift and VERITAS observations in addition to radio and optical monitoring 
by OVRO, MOJAVE, and CRTS. A strong flare was observed in optical, UV, and X- 
rays on 2012 December 30, quasi-simultaneously to the 7-ray flare, reaching a record 
flux for this source from optical to 7 rays. VERITAS observations at very high energy 
(E > 100 GeV) during 2013 January 6-17 resulted in an upper limit of E>o.2TeV < 
4.0 x 10~ 12 phcm -2 s _1 . We compared the spectral energy distribution (SED) of the 
flaring state in 2013 January with that of an intermediate state observed in 2011. The 
two SEDs, modelled as synchrotron emission and an external Compton scattering 
of seed photons from a dust torus, can be modelled by changing both the electron 
distribution parameters and the magnetic field. 

Key words: galaxies: active - galaxies: nuclei - galaxies: Seyfert - galaxies: individ- 
ual: PMN, J0948+0022 - gamma-rays: general 


1 INTRODUCTION 
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Relativistic jets are the most powerful manifestations of 
the release of energy produced by super-massive black holes 
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(SMBH) in the centers of active galactic nuclei (AGN). In 
about 15 per cent of AGN the accretion disc is at the base of 
a bipolar outflow of relativistic plasma, which may extend 
well beyond the host galaxy, forming the spectacular lobes of 
plasma visible in radio. The jet emission is observed across 
the entire electromagnetic spectrum. When the jet axis is 
closely aligned with our line of sight, the rest-frame radiation 
is strongly amplified due to Doppler boosting with a large 
fraction of the output observed at high energy, giving rise to 
the blazar phenomenon. Radio galaxies are typically viewed 
at larger angles, with a less ext reme amplification of the 
emission lilJrrv fc Padovaniill995 :i. As a consequence, a very 
small fraction of r adio galaxies is observed at high energy 
dAckermann et ahll201ll ). The discovery in 2008 by the Large 
Area Telescope (LAT) on board the Fermi satellite of 7-ray 
emission from radio-loud narrow-line Seyfert Is (NLSyl) re- 
vealed the presence of a possible new class of AGNs with 
relativistic jets. PMN J0948+0022 was the first NLSyl de- 
tected by Fermi - LAT dAbdo et al.ll2009a|). followed bv the 
detection of four other N LSyls in 7 rays l Abdo et al.i r2009c; 
ID’Ammando et alJl20l3 l. 

NLSyl is a class of AGN identified by 

lOsterbrock fe Poggel (1 19851 ) and defined in terms of 
the optical properties of its members: narrow permitted 
lines (FWHM (H/3) < 2000 k m s" 1 ) [OIII]/H/3 < 3, and 
a bump due to Fe II (see e.g. IPoggg 120001 . for a review). 
Variability and spectral properties from radio to 7-ray 
bands, together with the spectral energy distribution (SED) 
modeling, indicate a blazar-like behaviour of these indi- 
vidual objects (e.g. Abdoet al.l l2009cd : ID’Ammando et~ahl 
l2013al : IFoschini etal.1 20121 1. One of the most surprising 
aspects related to these NLSyls was the detection of flaring 
activity in 7 rays. A few st rong 7-ray flar e s were observed 
from PMN J0948+0022 dDonatol mid : iFoschinil 1201(1 
ID’Ammando fe Ciprinil 1201 lh . with an apparent isotropic 
7-ray luminosity up to 10 48 erg s _1 , comparable to that 
of the bright fl at spectr um radio quasars (FSRQ; e.g., 
PKS 1510—089; lOrienti et al.1 l2013l l. Similar 7-ray flares 
have been observed from two other NLSyls: SBS 0846+513 
( Donatol[201lj: U’A nnnando et al.ll2013al l. and 1H 0323+342 


( Carpenter et alF 20131 ). In addition, the high apparent 


7-ray luminosity observed in these objects suggests a large 
Doppler boosting and thus a small viewing angle, similar to 
blazars. These are important indications that NLSyls are 
able to host relativistic jets as powerful as those in blazars. 

PMN J0948+0022 could be considered the archetypal 
object of a new class of 7-ray emitting AGNs; therefore 
several multifrequency campaigns were performed to inves- 
tigate in detail its characteristics over the whole electro- 
magnetic spectrum ( ALdoetal. 2009b: IFoschini et al]l201lL 
Soil; [d: Ammando et al. 2014). New 7-ray flaring activity 
from PMN J0948+0022 was detec ted by Ferrm-LAT in 2013 
January (D’Ammando & Orientil 12013d ). triggering follow- 
up observations in optical, UV, and X-rays by Swift, and for 
the first time at very high energy (VHE; E > 100 GeV) by 
VERITAS. These data are complemented by the monitoring 
campaigns performed by the Owens Valley Radio Observa- 
tory (OVRO) and the Monitoring Of Jets in Active galactic 
nuclei with VLBA Experiments (MOJAVE0) at 15 GHz, 


1 http : / /www. physics .purdue . edu/MOJAVE/ 


and by the Catalina Real-time Transient Survey (CRT50) 
in V-band. 

The aim of the paper is to discuss the radio-to-7-ray ac- 
tivity from PMN J0948+0022 during 2012 December-2013 
August. This allows us to characterize the SED of the 2013 
flaring state of the source in order to compare the SED with 
that observed in a different activity state observed in 2011 
May, to investigate the emission mechanisms of this source 
and to discuss its properties in the context of the blazar 
scenario. 

The paper is organized as follows. In Sections 2 and 3, 
we report the LAT and VERITAS data analysis and results, 
respectively. The X-ray, UV, and optical data collected by 
Swift are presented in Section 4, together with the optical 
data collected by CRTS. Radio data collected by the OVRO 
40 m telescope and MOJAVE are presented in Section 5. In 
Section 6, we discuss the properties of the source and the 
modelling of the SEDs, and draw our conclusions in Section 
7. 

Throughout the paper the quoted uncertainties are 
given at the la level, unless otherwise stated, and the pho- 
ton indices are parameterized as dN/dE oc F -1 ’" , where r„ 
is the photon index at the different energy bands. We adopt 
a A cold dark matter (A-CDM) cosmology with H 0 = 71 km 
s _1 Mpc~ 4 , SJa = 0.73, and S2 m = 0 .27. The corresp onding 
luminosity distance at z = 0.5846 ((Schneider et al.l 12010 ). 
the source redshift, is dz, = 3413 Mpc, and 1 arcsec corre- 
sponds to a projected size of 6.59 kpc. 


2 FERMI - LAT DATA: ANALYSIS AND 

RESULTS 

The Fermi-LAT is a pair-conversion telescope operating 
from 20 MeV to > 300 GeV. It has a large peak effective 
area (~ 8000 cm 2 for 1 GeV photons), an energy resolution 
of typically ~10 per cent, and a field of view of about 2.4 sr 
with single-photon angular resolution (68 per cent contain- 
ment radius) of 0?6 at E = 1 Ge V on-axis. Further de tails 
about the Fermi-LAT are given in I Atwood et al.l d 20091 ). 

The LAT data reported in this paper were collected 
from 2012 December 1 (MJD 56262) to 2013 August 31 
(MJD 56535). During this time, the Fermi observatory op- 
erated almost entirely in survey mode. The analysis was 
performed with the ScienceTools software package ver- 
sion v9r32p5 0. Only events belonging to the ‘Source’ class 
were uscqj. The time intervals when the rocking angle of 
the LAT was greater than 52° were rejected. In addition, 
a cut on the zenith angle (< 100°) was applied to re- 
duce contamination from the Earth limb 7 rays, which are 
produced by cosmic rays interacting with the upper at- 
mosphere. The spectral analysis was performed with the 
instrument response functions P7REP_S0URCE_V15 using an 
unbinned maximum-likelihood method implemented in the 
tool gtlike. Isotropic (iso_source_v05.txt) and Galactic dif- 
fuse emission (gll_iem_v05_revl.fit) components were used to 


2 http://crts.caltech.edu/ 

3 http : //f erml . gsfc.nasa.gov/ssc/data/analysis/software/ 

4 http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/documentation/ 
Cicerone/Cicerone_L AT _IRF s /IRF_overview. html 
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model the background The normalizations of both compo- 
nents were allowed to vary freely during the spectral fitting. 

We analysed a region of interest of 10° radius cen- 
tred at the location of PMN J0948+0022. We evaluated the 
significance of the 7-ray signal from the source by means 
of the maximum-likelihood test statistic TS = 2 (logLi - 
logLo), where L is the likelihood of the data given the model 
with (Li) or without (Lp ) a point source at the position 
of PMN J0948+0022 (e.s.. iMattox et ahlll996l 'l. The source 
model used in gtlike includes all of the point sour ces from 
the second Fermi-LAT catalogue l2FGL: fNolan et al .11 2 01 j ) 
and the preliminary list of the third Fermi- LAT catalogue 
(3FGL; Ackermann et al., in prep.) that fall within 15° of 
the source. The spectra of these sources were parametrized 
by power-law functions, except for 2FGL J0909. 1+0121 and 
2FGL J1023. 6+0040 for which we used a log-parabola as in 
the 2FGL catalogue. A first maximum-likelihood analysis 
was performed in the 0.1-100 GeV energy range to find and 
remove from the model those sources that have TS < 10 
and/or the predicted number of counts based on the fitted 
model Npred < 3. A second maximum-likelihood analysis 
was performed on the updated source model. In the fitting 
procedure, the normalization factors and the photon indices 
of the sources lying within 10° of PMN J0948+0022 were 
left as free parameters. For the sources located between 10° 
and 15° from our target, we kept the normalization and the 
photon index fixed to the values from the 2FGL catalogue 
or from the 3FGL catalogue for those sources not reported 
in the former catalogue. 

Integrating over the period 2012 December 1 - 2013 
August 31 (MJD 56262-56535) in the 0.1-100 GeV 
energy range, using a log-parabola (LP), dN/dE oc 
(E/Eo)- a ~ p lo s (B/B o), as in the 2FGL catalogue, the fit 
yielded for PMN J0948+0022 a TS = 638, with an average 
flux of (15.5 ± 1.1) xlO -8 ph cm -2 s -1 , a spectral slope of a 
= 2.55 ± 0.08 at the reference energy Eo = 271 MeV, and a 
curvature parameter /3 = 0.05 ± 0.02. Using a simple power 
law (PL) we obtained a TS = 630, with a photon index of T 
= 2.62 ± 0.06 and an average flux of (15.9 ± 1.0) xl0~ 8 ph 
cm~ 2 s -1 . We used a likelihood ratio test (LRT) to check 
a PL model (null hypothesis) against a LP model (alterna- 
tive hypothesis). T hese values may be compared, following 
iNoIan et alJ (120121') . by defining the curvature test statis- 
tic TScurve = (TSlp - TSpl). The LRT results in a TS CU rve 
= 8, corresponding to a ~2.8 a difference; therefore a hint 
of curvature was observed in the 7-ray spectrum of PMN 
J0948+0022 during 2012 December-2013 August. By con- 
sidering only the period 2012 December 1 - 2013 June 14, 
in which the source was regularly detected on weekly time- 
scales, the maximum-likelihood analysis results for the LP 
and PL models in a TSlp = 677 and TSpl = 660, respec- 
tively. The corresponding curvature test statistic is TS CU rve= 
17, which corresponds to a significance of ~4.3cr. This indi- 
cates that the LP is preferred over a simple PL for this 
time interval, in agreement with the results obt ained for the 
obser vations during 2008 August-2010 August dNolan et ali 
120121+ 

Figure [I] shows the 7-ray light curve for the period 2012 
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Figure 1. Integrated flux light curve of PMN J0948+0022 ob- 
tained by Ferrai-LAT in the 0.1 100 GeV energy range during 
2012 December 1-2013 January 31 with 7-day time bins. Arrow 
refers to 2a upper limits on the source flux. Upper limits are com- 
puted when TS < 10. Open circles represent daily fluxes reported 
for the periods of high activity. No upper limits are shown for daily 
data. The horizontal line indicates the period of the VERITAS 
observation. 

December 1-2013 January 31 using a LP model and 1-week 
time bins. For the highest significance periods we also re- 
ported fluxes in 1-day time intervals. For each time bin, the 
spectral parameters for PMN J0948+0022 and for all the 
sources within 10° from the target were frozen to the value 
resulting from the likelihood analysis over the whole period. 
If TS < 10, 2 <t upper limits were calculated. All uncertain- 
ties in measured 7-ray flux reported throughout this paper 
are statistical only. The systematic uncertainty on the effec- 
tive area amounts to 10 per cent below 100 MeV, decreasing 
linearly with the logarithm of energy to 5 per cent between 
316 MeV and 10 GeV, and increasing line arly with the loga- 
rithm of energy up to 15 per cent at 1 TeV dAckermann et all 

l2012h . 

During the period of high activity (2012 December 15 
- 2013 January 11; MJD 56276-56303), the fit with a LP 
model in the 0.1-100 GeV energy range results in a TS = 
410, with an average flux of (46.8 ± 4.7) x 10~ 8 ph cm -2 s _1 , 
a spectral slope of a = 2.37 ± 0.10 at the reference energy 
Eo = 271 MeV, and a curvature parameter (5 = 0.10 ± 
0.04. This suggests marginal spectral variability during the 
7-ray flaring activity with respect to the whole period. The 
emission peak was observed on 2013 January 1 (MJD 56293), 
with an average flux for that day of (155 ± 31)xl0~ 8 ph 
cnU 2 s” 1 in the 0.1-100 GeV energy range, corresponding to 
an apparent isotropic 7-ray luminosity of ~1.5x 10 48 erg s _1 . 
This is the highest flux observed from PMN J0948+0022 so 
far in 7 rays (see iFoschini et, aTI 1 2 01 2l + 

By means of the gtsrcprob tool and using the source 
model optimized over the period 2012 December - 2013 Au- 
gust, we estimated that the highest energy photon detected 
from PMN J0948+0022 was observed on 2013 January 5 at a 
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distance of 0?06 from the source with a probability of 98.5% 
to be associated with the target and an energy of 7.2 GeV. 
Analyzing the LAT data collected over the whole period in 
the 1-100 GeV energy range with a simple power law, the 
fit yielded a TS = 143, with a photon index of F 7 = 2.64 
± 0.19 and an average flux of (3.4 ± 0.5) xl0~ 9 ph cm -2 
s _1 . On the contrary, the source is not detected at E > 10 
GeV (TS=1) with a 2cr upper limit of 1.5xlO -10 ph cm -2 
s' 1 (assuming a photon index T 7 =3). 


3 VERITAS DATA: ANALYSIS AND RESULTS 

The VERITAS observatory Jllolder et al .1 [200(3 . l2008h is an 
array of four imaging atmospheric Cherenkov telescopes 
(IACT) located at the Fred Lawrence Whipple Observa- 
tory near Tucson, Arizona. Each telescope consists of a 12- 
m diameter reflector and a photomultiplier camera covering 
a field of view of 3? 5. The array has an effective area of 
~5x 10 4 m 2 between 0.2 and lOTeV. 

VERITAS can detect a source with 2-3 per cent of the 
Crab Nebula fltu0 in 300 min of exposure, with improved 
sensitivity to soft-spectrum sources following an up date of 
the c amera photomultipliers in Summer 2012 llKieda et al.l 
l2013l !. The angular and energy resolution for reconstructed 
7-ray showers are 0? 1 and 15 per cent, respectively, at 1 TeV. 

After a tr igger from Fermi- LAT 

dD’Ammando fc Orient! l2013ch . VERITAS observations 
of PMN J0948+0022 took place between 2013 January 6 
(MJD 56298) and January 17 (MJD 56309), covering the 
decay of the 7-ray flare. Observations consisted of 30 min 
exposures pointing 0? 5, offset in wobble observing mode 
dFomin et al.lll994 l at a median zenith angle of 32°. About 
one third of the collected data had to be rejected because 
of cloud coverage, leaving an effective dead-time corrected 
exposure of 315 min. 

A description o f the VERITAS analysis can be found in 
lAcciari et al J d2008l f and I Archambault et aid (120131 ) . To en- 
sure a robust estimation of the shower parameters, at least 
two telescope images with a total signal size equivalent to 
> 74 photoelectrons were required for an event to be recon- 
structed. Background-rejection cuts, optimized for a mod- 
erately bright soft-spectrum source, were applied to remove 
over 99.9 per cent of the charged cosmic rays. Events with 
reconstructed arrival direction 9 < 0?1 to the position of 
PMN J0948+0022 were selected as signal. The number of 
background events in the signal region was estimated from 
off-source regions in the same field o f view using the reflected 
region technique llBerge et al.ll2007h . 

The analysis revealed no significant 7-1'ay signal from 
PMN J0948+0022. An excess (A ex ) of 11 candidate 7-ray 
events over a background of 615 was measured, correspond- 
ing to a si gnificance of 0.4 standard deviations according to 
Eq. 17 in lLi fc Mai (jl983l L with N on = 626, N 0 g = 5211, 
and a ratio of the on-source time to the off-source time 
a = 0.118. Assuming a power-law spectrum with photon 
index T 7 = 3.5, a 99 per cent confidence level (c.l.) flux up- 
per limit of E>o.2TeV < 4.0 x 10 -12 phcnD 2 s _1 was derived 


6 1 Crab = 2.1 X 10~ 10 

llHillas et alj|T993l 


ph c 


using the prescription bv lRolke et al.l J2005l l. The systematic 
uncertainty in the flux measurement is approximately 40 per 
cent. Differential upper limits (95 per cent c.l., T 7 = 3.5 to 
match the LAT spectrum during the period of high activ- 
ity) at the energy threshold of the analysis (170 GeV) and 
where the correlation between flux level and spectral shape 
is minimal (300 GeV) are shown in Fig. [3] Compatible re- 
sults were obtained using an independent analysis package 
with a different set of event-selection cuts. 

Given the variable nature of the 7-ray emission from 
PMN J0948+0022, a search for signal in nightly bins (Ta- 
blc[TJ as well as in 30-min runs was also conducted. However, 
no significant 7-ray excess was found at those time-scales ei- 
ther. 


4 X-RAY TO OPTICAL DATA 


4.1 Swift data analysis 

The Swift satellite (iGehrels et al.ll2004f) performed six obser- 
vations of PMN J0948+0022 between 2012 December and 
2013 January. The observations were carried out with all 
three instruments o n board: the X-ray Telescope (XRT; 
iBurrows et al.ll2005l . 0.2-10.0 ke V), th e Ultraviolet /Optical 
Telescope fUVOT: lRoming et al.ll2 005l. 170-600 nml and the 
Burst Alert Telescope ('BAT; iBarthelmv et al.ll2005l . 15-150 
keV). The hard X-ray flux of this source is below the sensi- 
tivity of the BAT instrument for the short exposures of these 
observations, therefore the data from this instrument are not 
used. Moreover, the source was n ot present in the Swift BA T 
70- month hard X-ray catalogue dBaumgartner et al.l 2013j j . 

The XRT data were processed with standard procedures 
(xrtpipeline vO.12.8), filtering, and screening criteria us- 
ing the HEAsoft package (v6.13). The data were collected 
in photon counting mode for all of the observations. The 
source count rate was low (< 0.5 counts s^ 1 ); thus pile-up 
correction was not required. Source events were extracted 
from a circular region with a radius of 20 pixels (1 pixel = 
2.36 arcsec), while background events were extracted from a 
circular region with radius of 50 pixels away from the source 
region and from other bright sources. Ancillary response files 
were generated with xrtmkarf , and account for different ex- 
traction regions, vignetting and point-spread function cor- 
rections. We used the spectral redistribution matrices v013 
in the Calibration data base maintained by HEASARC0. 
The spectra with low numbers of photons collected (< 200 
counts) were rebinned with a min imum of 1 count per bin 
and the Cash statistic (|Cas IB was used. We fitted the 
spectra with an absorbed power-law using th e photoelectric 
absorption model tbabs (I Wilms et alj 120001 ) . with a neu- 
tral hydrogen co lumn density f i xed t o its Galactic value 
(5.07 x 10 2 °cm~ 2 ; iKalberla et al.ll2005l l. The fit results are 
reported in Tabic El and the 0.3-10 keV fluxes are shown in 

Fig. El 

UVOT data in the v, b, u, wl, m2, and w2 filters were 
reduced with the HEAsoft package v6.13 and the 20120606 
CALDB-UVOTA release. We extracted the source counts 
from a circle with 5 arcsec radius centred on the source and 
the background counts from a circle with 10 arcsec radius in 


for E > 200 GeV 
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Table 1. VERITAS observations and results of PMN J0948+0022. 


Date 

(MJD) 

Date 

CUT) 

Exposure 

(min) 

Significance 

(Vex 

F> 0.2TeV 

(10 -11 phcm -2 s -1 ) 

56298 

2013-01-06 

33 

-1.2 

-9.8 

< 0.84 

56299 

2013-01-07 

77 

-0.3 

-3.3 

< 0.54 

56301 

2013-01-09 

75 

1.1 

15.0 

< 1.12 

56305 

2013-01-13 

25 

0.0 

-0.1 

< 1.22 

56306 

2013-01-14 

17 

0.0 

0.0 

< 1.53 

56307 

2013-01-15 

13 

1.9 

10.6 

< 4.58 

56308 

2013-01-16 

50 

0.3 

3.0 

< 1.08 

56309 

2013-01-17 

25 

-1.2 

-8.8 

< 2.22 

all data 

- 

315 

0.4 

10.6 

< 0.40 


Table 2. Log and fitting results of S‘imJi-X\{T observations of PMN J0948+0022 using a power-law model with a HI column density fixed 
to the Galactic value in the direction of the source. “Unabsorbed flux. 


Date 

(MJD) 

Date 

(UT) 

Net exposure time 
(8) 

Photon index 
(Fx) 

Flux 0.3-10 keV“ 
(10 -12 erg cm -2 s -1 ) 

56264 

2012-12-03 

1963 

1.53 ± 0.25 

3.3 ±0.5 

56286 

2012-12-25 

1893 

1.51 ± 0.21 

5.1 ±0.5 

56291 

2012-12-30 

722 

1.77 ± 0.16 

13.6 ±1.3 

56295 

2013-01-03 

2972 

1.48 ± 0.15 

6.8 ±0.5 

56303 

2013-01-11 

2932 

1.60 ± 0.17 

6.9 ±0.6 

56309 

2013-01-17 

3284 

1.67 ± 0.10 

10.6 ±0.6 


a nearby source- free region. The o bserve d mag nitudes are 
reported in Table [3] As in, e.g., iRaiteri et al.l (1201 lh . we 
calculated the effective wavelengths, count-to-flux conver- 
sion factors (CFa), and amount of Galactic extinction in 
the UVOT ban ds (Aa) by con volving the mean Galactic ex- 
tinction law bv ICardelli et all (| 19891 ) with a power-law fit to 
the source flux and with the filter effective areas. To obtain 
the UVOT de-reddened fluxes we multiplied the count rates 
for the CFa and corrected for the corresponding Galactic ex- 
tinction values Aa- The UVOT density fluxes are reported 
in Fig H 


4.2 CRTS data 

The source has been monito red by the CRTS dDrake et al.l 
120091 : iDiorgovski et al.l|201ll ) , using the 0.68-m Schmidt tele- 
scope at Catalina Station, AZ, and an unfiltered CCD. The 
typical cadence is four exposures separated by 10 min in a 
given night; this may be repeated up to four times per lu- 
nation, over a period of ~6-7 months each year, while the 
field is observable. Photometry is obtained using the stan - 
dard Source- Extractor package (IBertin fe Arnouta 1 1 99(3) . 
and transformed from the unfiltered instrumental magni- 
tude to Cousins V0 by V = Vcss + 0.31 [B — V ) 2 + 0.04. 
We averaged the values obtained during the same observing 
night. The flux densities collected by CRTS in V band are 
reported in Fig. [2] 


5 RADIO DATA 

5.1 OVRO data analysis 

As part of an ongoing blazar monitoring program, the 
OVRO 40-m radio telescope has observed PMN J0948+ 0022 
at 15 GHz regularly since the end of 2007 iRichards et ali 
Im3). This monitoring program includes about 1700 known 
and likely 7-ray loud blazars above declination —20°. The 
sources in this programme are observed in total intensity 
twice per week with a 4 mjy (minimum) and 3 per cent 
(typical) uncertainty on the flux densities. Observations 
were performed with a dual-beam (each 2.5 arcmin FWHM) 
Dicke-switched system using cold sky in the off-source beam 
as the reference. Additionally, the source is switched be- 
tween beams to reduce atmospheric variations. The absolute 
flux density scale is calibrated using obs ervations of 3 C 286 , 
adopting the flux density (3.44 Jy) from lBaars et all d 19771 ). 
This results in about a 5 per cent absolute scale uncertainty, 
which is not reflected in the plotted errors. 

5.2 MOJAVE data analysis 

We investigated the parsec-scale morphology and flux den- 
sity variability at 15 GHz by means of 6-epoch Very 
Long Base line Array data from the MOJAVE programme 
dLister et al.ll2009l ). The data sets span the time interval be- 
tween 2012 July and 2013 July, in order to overlap with the 
Fermi - LAT data. We imported the calibrated uv data into 
the National Radio Astronomy Observatory AIPS package. 
In addition to the total intensity images, we produced the 
Stokes Q and U images, to derive information on the po- 


http://nesssi.cacr.caltech.edU/DataRelease/FAQ2.html#improvelarized emission. The flux density was derived by means of 
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Table 3. Results of the Swift - UVOT observations of PMN J0948+0022 in magnitudes. 


Date (MJD) 

Date (UT) 

V 

b 

u 

wl 

m2 

w2 

56264 

2012-12-03 

- 

- 

17.43i0.03 

- 

- 

- 

56286 

2012-12-25 

18.lli0.18 

18.29i0.10 

17.34i0.07 

17.29i0.07 

17.46i0.05 

17.46i0.04 

56291 

2012-12-30 

16. OliO. 10 

16.59i0.07 

15.86i0.06 

16.20i0.04 

15.88i0.08 

16.08i0.05 

56295 

2013-01-03 

17.48i0.09 

18.00i0.07 

17.19i0.06 

17.07i0.05 

17.16i0.06 

17.12i0.04 

56303 

2013-01-11 

17.70i0.ll 

17.85i0.06 

17.02i0.05 

16.94i0.05 

17.01i0.05 

17.14i0.04 

56309 

2013-01-17 

17.45i0.09 

17.81i0.06 

16.99i0.05 

17.00i0.05 

17.02i0.05 

17.04i0.04 



56300 56400 56500 

Time [MJD] 


Figure 2. Multifrequency light curve for PMN J0948+0022. The period covered is 2012 December 1 - 2013 August 31 (MJD 56262- 
56535). The data were collected (from top to bottom) by Fermi-LAT (7 rays; in units of 10 -8 ph cm -2 s _1 ), Swift-XRT (0.3-10 keV; 
in units of 10 “ 12 erg cm -2 s~ 1 ), Swift - UVOT (rul. m2, and w2 filters; in units of mjy), CRTS (open circles; in units of mjy) and 
Swift - UVOT (?;, b, and u filters; in units of mjy), MOJAVE (open circles) and OVRO (15 GHz; in units of Jy). The horizontal line in 
the top panel indicates the period of the VERITAS observation. 


the AIPS task JMFIT which performs a Gaussian fit on the 
image plane. Total intensity flux density and polarization 
information are reported in Table [4] The uncertainties on 
the flux density scale are less than 5 per cent, a nd the errors 
on the polarization angle are ^5° IILister et al.ll2013i) . 


6 DISCUSSION 

6.1 Radio to 7-ray behaviour 

During 2012 December - 2013 January the NLSyl PMN 
J0948+0022 showed new flaring activity in 7 rays, triggering 
follow-up observations by Swift and VERITAS. In Fig. [2] we 
compare the 7-ray light curve collected by Fermi - LAT in the 
0.1 TOO GeV energy range with the X-ray (0.3-10 keV), UV 
(wl, m2, and w2 filters), optical ( v , b, and u filters), and 
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Table 4. Flux density and polarization of PMN J0948+0022 from 
15 GHz MOJAVE data. 


Date 

(MJD) 

Date 

(UT) 

Score' 1 

( m Jy) 

Sjet a 

(mjy) 

S po i a 

(mjy) (%) 

EVPA 6 

(deg) 

56120 

2012-07-12 

325 

5 

3 

(0.9%) 

96 

56242 

2012-11-11 

336 

5 

3 

(0.9%) 

86 

56271 

2012-12-10 

380 

6 

4 

(1.0%) 

97 

56313 

2013-01-21 

483 

4 

3 

(0.6%) 

130 

56417 

2013-05-05 

862 

5 

6 

(0.7%) 

96 

56503 

2013-07-30 

420 

3 

12 

(2.9%) 

90 


“Uncertainties on the flux densities are less than 5 per cent. 
^Uncertainties on the polarization angles are 5°. 


radio (15 GHz) light curves collected by Swift (XRT and 
UVOT), CRTS, MOJAVE, and OVRO. 

The 7 -ray flaring period peaked on 2013 January 1 
(MJD 56293) with a variability amplitude (calculated as the 
ratio between the maximum and minimum flux) of ~15 be- 
tween 2012 December and 2013 August, and flux variations 
on a daily time-scale. The peak flux reported here, (155 ± 
31)xl0 -8 ph cm -2 s , is the g reatest detected from this 
source so far at 7 rays (see e.g.. iFoschini et all [20121. The 
7 -ray flaring activity ended in 2013 mid-January, and no fur- 
ther significant increase of the 7 -ray flux was observed up 
to 2013 August. Only marginal spectral variability was ob- 
served d uring this flarin g activ ity, similar to the behaviour 
shown in IFoschini et al] (|2012| ) for the period 2008 August 
- 2011 December. 

During the highest 7 -ray activity a quasi-simultaneous 
increase was observed from optical to X-rays. Swift /X RT ob- 
served the source with a 0.3-10 keV flux in the range (3.3- 
13.6) xl0 ~ 12 erg cm ~ 2 s -1 , with a variability amplitude of 

(Table 0 • The peak of the X-ray flux was observed on 
2012 December 30 (MJD 56291), when PMN J0948+0022 
reached the highes t flux detected for this source so far (for 
a com parison, see IFoschini et al.l 120121 : iD’Ammando et~ahl 
1201411 . No significant spectral change in the X-ray band was 
observed during the Swift observations, with the photon in- 
dex ranging bet ween 1.5 and 1.8, similar ly to what was ob- 
served in 2011 (ID’Ammando et al.ll2014h . It is worth men- 
tioning that a second increase of the X-ray flux was observed 
on 2013 January 17 (MJD 56309), but no corresponding in- 
crease of activity was observed in 7 -ray, UV, and optical 
bands. 

The optical and UV emission significantly increased 
over only 5 days, from 2012 December 25 (MJD 56286) 
to December 30 (MJD 56291), suggesting that the domi- 
nant contribution to the continuum flux should come from 
the synchrotron emission. However, the variability ampli- 
tude decreases with frequency: from ~7 in V-band to ~3.5 
in the UV part of the spectrum. This could be due to the 
contribution in UV of the thermal emission from the accre- 
tion disc, whi ch is well identified d uring low activity periods 
of this source (IFoschini et al.|[2012f) . A similar behavi our was 
alrea dy noticed in some FSRQs (e.g. 3C 454.3. IRaiteri et alJ 
EH' After the 7 -ray flare, two observations were performed 


in V-band by the Steward Observatory as part of the blazar 
monitoring programme of the University of Arizon43 Optical 
polarization of 4.0±0.1 per cent and 1.1±0.1 per cent was 
estimated on 2013 January 17 and 18, respectively, while the 
position angle on the sky of the linear polarization changed 
from 100°± 1° to 26°± 3° . The value observed on 2013 Jan- 
uary 17 is not as high a s the 12 .3 per cent observed for this 
source on 2013 May 25 dEggen et al.ll20l3 l. but the change 
of a factor of ~4 in one day is notable as well as the signif- 
icant change of the observed polarization angle. The rapid 
variation of optical polarization is a characteristic associ- 
ated with the variation of synchrotron emission in blazars. 
The variation of both the fractional polarization and polar- 
ization angl e may b e an indication of a turbulent magnetic 
field fe.g.. lMarscheill2014l ). 

PMN J0948+0022 was highly variable at 15 GHz dur- 
ing the OVRO 40 m telescope monitoring, with the flux 
density varying from 357 mjy (on MJD 56281; 2012 De- 
cember 20) to 936 mjy (on MJD 56445; 2013 June 2), as 
shown in Fig. [2] The source was observed in a low state 
at 15 GHz during the peak of the 7 -ray activity. On the 
other hand, no significant 7 -ray activity was detected at the 
time of the peak of the radio activity. If the 7 -ray and ra- 
dio activity are connected, the observed peak at 15 GHz is 
produced by the same disturbance of the jet responsible for 
the 7 -ray flare delayed by ~5 months. This delay may be 
explained by synchrotron self-absorption opacity effects in 
combination with shock propagation effects. As an alterna- 
tive, the 7 - 1 'ay and radio emission may originate in different 
parts of a bending and inhomogeneous jet, with variable 
orientations with respect to the line of sight. Therefore we 
may have different Doppler boostings of the zone respon- 
sible for the optical-to- 7 -ray emission and that responsi- 
ble for the radio e mission in the two flaring periods (e.g., 
IRaiteri et alJl20l3 ). A delay of ~50 days of the 15 GHz flux 
density peak with respect to the 7 -ray one was observed 
during the 2010 7 -ra y flaring activity of PMN J0948+0022 
(IFoschini et al.l l2011). On the contrary, no obvious connec- 
tion between 7 and radi o emission was observed in 2011 
(jD’Ammando et al. I EH). The lack of a coherent pattern 
suggests a complex connection between the radio and 7 -ray 
emission for PMN J0948 +0022, as also observed in the ot her 
NLSyl SBS 0846+513 dD’Ammando et a l.ll2013a l. 1201211 a s 
well as in blazars (e.g. PKS 1510—089, Orienti et. al .1120 Ijl ). 

At pc-scales the radio emission is dominated by the core 
component, which accounts for ~98 per cent of the total 
emission. A hint of the jet emerges from the core with a 
position angle ^30 ° , consistent with what was found in pre- 
vious works (e.g., ID’Ammando et ~ahl 1201 ll : iGiroletti et al.l 
EH)- From the pc-scale resolution images we can separate 
the core flux density from the jet emission, in order to study 
the light curve and polarization properties of each of them. 
In accordance with the OVRO monitoring, the flux density 
of the core observed by MOJAVE at 15 GHz increased from 
2012 July 12 to 2013 May 5, when a value of 862 mjy was 
reached (Table|4|. No significant increase of the flux density 
from the jet was observed. A higher polarized emission of 
the core (S'pol) and polarization percentage were observed 
on 2013 July 30, after the peak of the radio emission. The 


9 


http : / /james . as . arizona.edu/~psniith/Feriiii/ 


8 F. D’Ammando, M. Orienti, J. Finke, et al. 


electric vector position angle (EVPA) of the core shows a 
moderate change during the MOJAVE observation, ranging 
from 86 ° to 130 °, similar to what was ob served for the same 
object in 2011 (ID’Ammando et al .11201 -'ll ). In particular, we 
noted a change of ~25° between 2012 December and 2013 
January, after the 7 -ray flare observed by LAT. 


Because 7 -ray emission from NLSyls was not expected 
until its discovery with the Fermi satellite, no dedicated 
observations of these sources were carried out at VHE be- 
fore that time. The detection of a few of these objects in 7 
rays by the LAT has given rise to a new interest in their 
emission processes at the higher energies. The only 7 -ray 
emitting NLSyl with VHE observations in the literature is 
1H 0323+342 (z = 0.062), which was selected from a list 
of candidate TeV-emitting FSRQs (IPerlmani l2000l l . It was 
observed by the Whipple telescope during 2001-2003, re- 
sulting in an upper limit of 5.2xl0 ~ 12 ph cm -2 s~ x above 
400 GeV, with margina l evi d ence o f a ra te increase in the 
2001-2002 light curve (iFalcone et al .1 [20041) . This paper pro- 
vides the first report of an observation of a NLSyl with 
the current generation of IACTs (i.e. VERITAS, MAGIC, 
and H.E.S.S.). Unfortunately the source was not detected 
by VERITAS at VHE in 2013 January, with an upper limit 
of 4.0xl0 -12 ph cm -2 s _1 above 200 GeV. This could be 
due to several reasons. First of all, the distance of the source 
is relatively large ( z = 0.5846), and so most of the GeV /TeV 
emission should be absorbed due to pair production from 7 - 
ray photons of the source and the infrared photons from the 
extragalactic background light (EBL), although the mos t 
distant FSRQ detected at VHE, 3C 279 ll Albert et al.ll2008lh 
is at a comparab le distance of z = 0.5362. Considering the 
1 Finke et all (i2010l l model for the EBL the optical depth is r 
~1 for 170 GeV photons at z = 0.5846. As a comparison, us- 
ing the lDominguez et al j (2011) model the optical depth is r 
~1 for 200 GeV photons at the same distance. However, the 
VERITAS observations were carried out a few days after the 
peak of the 7 -ray activity from PMN J0948+0022 detected 
by LAT, thus covering only the last part of the MeV-GeV 
flare. An extrapolation to the VHE band of the LAT mea- 
sured energy spectrum averaged over the active phase of 
the flare (2012 December 29 - 2013 January 18) falls be- 
low the VERITAS sensitivity in 5h of observation once the 
EBL absorption is taken into account (Fig. [3]). Assuming 
no additional spectral curvature than that already seen in 
the LAT energy range (LP with /3 = 0.1) and the flux sup- 
pression induced by the EBL, detectable VHE emission from 
PMN J0948+0022 would be expected during GeV flares with 
Fo. l—ioo GeV + 2xl0 ~ 6 ph cm -2 s _1 . Given the similarities 
between PMN J0948+0022 and 7 -ray detected FSRQs, the 
presence of a bright broad-line region (BLR) could produce 
a spectral break due to pair production, suppressing the flux 
beyond a few GeV and preventing a VHE detection. How- 
ever, the role of BLR photons in the emission/absorption of 
VHE 7 -rays is far from understoo d, and the de tections of the 
FSRQ 3C 279 (|Albert et all200slb 4G +21.35 (fAleksic et al.1 
l201ll L and PKS 1510-089 (Abramowski et ~alT 20131 1 have 
shown that the spectrum of some FSRQs extends to VHE 
energies during certain flaring periods. This may be an in- 
dication that during those high activity periods the 7 -ray 
emission was produced outside the BLR. 
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Figure 3. SED of PMN J0948+0022 in the MeV-to-TeV energy 
range. The Fermi-LAT spectrum was extrapolated to the TeV 
ene rg ies and corre cted for EBL absorption using the model of 
iFinke et al.1 1 2 0 1 fl l. Fermi- LAT and VERITAS data points and 
upper limits are shown. VERITAS upper limits are calculated at 
the energy threshold of the analysis and where the correlation 
between flux level and spectral shape is minimal. 


6.2 SED Modelling 


We created two SEDs for comparing the 2013 January 7 -ray 
flaring activity (presented in thi s paper) to the 2011 Ma y 
intermediate state (presented in ID’Ammando et alJl2014h . 
The 2013 flaring SED includes the LAT spectrum built 
with data centred on 2012 December 29-2013 January 18 
(MJD 56290-56310), optical, UV, and X-ray data collected 
by Swift (XRT and UVOT) on 2012 December 30 (MJD 
56291), and radio data collected by OVRO at 15 GHz on 
2013 January 3 (MJD 56295). The 2011 intermediate SED 
was built with LAT data centred on 2011 May 22- June 11 
(MJD 55703-55723), optical, UV, and X-ray data collected 
by XMM-Newton (EPIC and Optical Monitor) on 2011 May 
28-29 (MJD 55709-55710), and radio data collected by Ef- 
felsberg at 15 GHz and 32 GHz on 2011 May 24 (MJD 
55705). 

The location along the jet of the 7 -ray emitting region 
is controver sial; see the discussion by Finke] d2013b| b Rapid 
variability llTavecchio et a.1.1 120101 : iBrownl 120131 : ISaito et al.l 
l2013ll indicates a compact emitting region that, if the emit- 
ting region takes up the entire cross section of the jet, 
argues for emission from inside the BLR. Absorption fea- 
tures in 7 -ray spectra could also be an indicat i on o f this 
dPoutanen fc Sternl2010l ; [stern fe Poutanenl201ll . l2014jh On 
the other hand, 7 -ray spectra that exten d to V HEs with- 
out absorption features in some objects llAleksic et al]l20lH : 
lAleksic et al.l 120141 : IPacciani et al.l l2014h are a strong indi- 
cation that emission takes place outside the BLR, where the 
dust torus is likely the dominant external radiation held. 
This is corroborated by the association of 7 -ray hares with 
the radio outbursts and the ejection of superluminal com- 
ponents from cores of blazars as observed with very long 
baseline interferometry, giving distances > a few pc based 
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on light travel tim e arguments (e.g., iMarscher et al.ll20ld : 
IQrienti et al.ll2013l ). There is some indication that different 
flares, e ven from the same source, can occur in different lo- 
cations llBrownll2013l : iNalewaiko et al.ll2014h . Equipartition 
argu ments can also be in voked for locating the emitting re- 
gion dPermer et al.1 12014* '). In the following, we first model 
the emitting region consistent with a dust torus external 
radiation field, then discuss modelling consistent with an 
emitting region inside the BLR. 

We modelled the SED of the two activity states of 
PMN J0948+0022 with a combination of synchrotron, syn- 
chrotron self-Compton (SSC), and external Compton (EC) 
non-thermal emission. The synchrotron component consid- 
ered is self-absorbed below 10 11 Hz and thus cannot re- 
produce the radio emission. This emission is likely from 
the su p erpos ition of multiple self-absorbed jet components 
IIKonigi l 198111 . We also included thermal emission by an ac- 
cretion disc and dust torus . Th e modelling de t ails c an be 
found in iFinke et al.l (120081 1 and IDermer et all (l2009h . Ad- 
ditionally, a soft excess was observed in the X-ray spectrum 
during the 2011 state, w hose origin is still under debate 
dD’Ammando et al.li20 14). In order to account for the soft 
X-ray excess in the 2011 SED, we included emission from 
the disc which is Compton scattered by an optically thin 
thermal plasma near the accretion disc (i.e., a corona). Thi s 
was done using the routine “SIMPL” (ISteiner et al.l 120091) . 
This routine has two free parameters: the fraction of disc 
photons scattered by the corona ( f sc ), and the power-law 
photon index of the scattered coronal emission (T sc ). The 
mass of the BH was c hosen to be the same as the one used by 
iFoschini et a ] d2Qll|), M bh = 1.5 x 10 8 Mm, cons istent with 
estimates by Zhou et "all d2003l l and lAbdo et al.l d2009al l. 

The observed variability time-scale of ~1 day (Fig. [T]) 
constrains the size of the emitting region, consistent with the 
7-ray light curve (Fig. [2]). The results of the modelling can be 
found in Table [5] mid Figure [4] (f or a description of the model 
parameters see [Dermer(et al. 2009). Both models have the 
jet power in magnetic field about a factor of 10 above the jet 
power in electronic. However, it is possible that the com- 
bined power in electrons and protons in the jet could be in 
equipartition with the magnetic field. If this is the case, it 
implies that the jet power in protons is Py p « 4.8 x 10 45 
erg s _1 from the 2011 model and Py p « 1.9 x 10 46 erg s _1 
from the 2013 flaring model. This model assumes the emit- 
ting region is outside the BLR, where dust torus photons 
are likely the seed photon source. This source was mod- 
eled as being an isotropic, monochromatic radiation source 
with dust parameters chosen to be consistent with the re- 
lation between inner radius, d i sc lum inosity, and dust tem- 
perature from [Nenkova et alj d2008h . Note that the model 
parameters shown here are not unique, and other model pa- 
rameters could also reproduce the SED, especially for the 
2011 state. 

In 2011 May- June the 7-ray and X-ray flux was inter- 
mediate between the 2013 flaring activity a nd the low ac- 
tivity state observed e.g. in 2009 May 15 dFoschini et al.1 
I2012f) . with no corresponding flare in the optical. This as- 
pect and the shape of the optical spectrum indicate that 


10 Jet pow e rs we re calculated assuming a two-sided jet. See 
iFinke et all DoOtj ). 


the optical emission during that period was dominated by 
thermal accretion disc emission. The synchrotron compo- 
nent was poorly constrained, other than an upper limit from 
the optical data. According to our model, thermal emission 
from the dust torus would be visible if infrared observations 
were available. The X-rays originated from a combination of 
coronal emission (for the soft excess) and EC emission. In 
this model, the SSC emission was quite weak and did not 
play a significant role in the SED. 

The 2013 activity state was brighter across the elec- 
tromagnetic spectrum than during 2011. In particular, the 
2013 flaring period included a flare in optical as well as 
in X-rays and 7 rays. The much stronger optical emission 
implies that the two activity states differed by more than 
just the electron distribution. This is in contrast to the 
FSRQ PKS 0537—441, where high and low activity states 
were modelled by only ch anging the electron distribution 
dD’Ammando et al.ll2013blh In this way, PMN J0948+0022 
is similar to the FSRQ PKS 2142—75 dDutka et al.ll2013h . 
which exhibited multiple flares where another parameter be- 
sides the electron distribution needed to be changed between 
flares. As for PKS 2142—75, we modelled the two activity 
states of PMN J0948+0022 by changing the magnetic field 
and the electron distribution. This may be an indication that 
the shock responsible for the flare produces a compression 
as it passes through the emitting region yielding an amplifi- 
cation of the magnetic field during the flare. The 2011 inter- 
mediate state from PMN J0948 +0022 is similar to “flare B” 
from PKS 2142—75 (iDutka et al.ll2013h . since in both cases 
the optical part of the spectrum is dominated by the ther- 
mal disc emission. The 2013 flare from PMN J0948+0022 is 
most like “flare A” from PKS 2142—75 dDutka et al.ll2013l ). 
since for both of these flares the optical spectrum of the 
two sources is dominated by non-thermal synchrotron emis- 
sion. We also modeled the 2013 flaring state with parameters 
consistent with Compton-scattering of BLR Ha (6563 A) 
line radiation; again, see Figure 4 and Table 5. The model 
reproduces the data approximately equally as well as the 
dust scattering model. Thus, it is not possible to distinguish 
between the two scenarios by snapshot spectral modeling. 
However, we note the BLR scattering model is farther from 
equipartition, so this may be a weak argument in favor of the 
dust scatte ring scenario. As also obs erved in the NLSyl SBS 
0846+513 dD’Ammando et al.ll2012lL PMN J0948+0022 ex- 
hibits behavior similar to other FSRQs. The X-ray emission 
is dominated by EC, and the coronal emission is not visible, 
so that there is no evidence for a soft excess from the data 
in 2013. However, this could be due to the low sensitivity of 
Swift-XRT at ~1 keV. 

A clear correlation between Compton dominance^]] and 
the rest-frame peak sy nchrotron frequency {i' p f) for blazars 
was reported in[Finkc d2013al h This correlation is related to 
the contribution of the EC component, which results to be 
higher for larger values of Compton dominance. In Fig. [5] the 
Compton dominance versus ufj, is shown for sources from the 
second LAT AGN Catalog 12LAC: lAckermann et all 1 201 il l. 
In addition we plot the low and high activity sta tes of the 
NLSyl SBS 0846+513 dD’Ammando et al.ll2013ah and the 


11 Compton dominance is the ratio of the peak Compton lumi- 
nosity to peak synchrotron luminosity. 
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Figure 4. SEDs and models for the 2013 and 2011 activity states 
from PMN J0948+0022. The filled circles are the data from the 
2013 flaring state, and the open squ ares are the data f r om th e 
2011 intermediate state taken from iD’Ammando et al.l d2014l) . 
The dashed curve indicates the disc and coronal emission, and 
the dotted line indicates the thermal dust emission. Solid lines 
represent models consistent with scattering dust torus radiation, 
while the dashed-dotted curve represents a model consistent with 
the scattering of BLR Ho; (6563 A) radiation. Arrows refer to 
2cr upper limits on the source flux. The VERITAS uppe r limits 
are co rrected for EBL absorption using the model of lFinke et al.l 

teoich . 



Figure 5. The Compton dominan ce versus peak synchrot ron fre- 
quency for sources from the 2LAC jAckermanr^e^al ■ 20111). along 
with two states of SBS 0846+513 jD’Ammando et al.ll2013ah and 
PMN J0948+0022 (this work). Filled circles represent FSRQs, 
empty circles represent BL Lac objects, an d filled square s repre- 
sent AGN of unknown type. Adapted from Finkel d2013al ) . PMN 
J0948+0022 and SBS 0846+513 are plotted as triangles. 


intermediate and high activity states of PMN J0948+0022 
shown in Fig. [4] Both NLSyls are in a region occupied by 
FSRQs, indicating the similarity between FSRQs and these 
NLSyls. It is worth mentioning that for both NLSyls the 
peak frequency is higher for lower Compton dominances, in 
agreement with the general behaviour observed for blazars. 


7 CONCLUSIONS 

We reported on the observation by the Fermi - LAT of flar- 
ing 7-ray activity from the NLSyl PMN J0948+0022 in 

2012 December-2013 January. The 7-ray flux reached on 

2013 January 1, (155 ± 31) xlO -8 ph cm -2 s -1 in the 0.1 
100 GeV energy range, is the highest value detected from 
this source and more generally from a NLSyl so far. Only 
marginal spectral variability was observed in 7 rays during 
the flare. PMN J0948+0022 was not detected above E > 10 
GeV by LAT. 

We presented multiwavelength observations of the source 
from radio to VHE during the period 2012 December 1-2013 
August 31 including data collected by VERITAS, Fermi- 
LAT, Swift, CRTS, MOJAVE, and OVRO. An intense flar- 
ing episode, quasi-simultaneous to the 7-ray one, was ob- 
served on 2012 December 30 in optical, UV, and X-rays, 
showing record flux also in these energy bands. This sug- 
gests a strict connection between the region responsible for 
the 7-ray activity and that responsible for the optical-to-X- 
ray activity, although the time sampling of the observations 
is too sparse from optical to X-rays to formally establish a 
correlation. Flaring activity was also observed at 15 GHz 
with the radio peak detected on 2013 June 2. The most 
likely scenario is that the radio flare is the delayed coun- 
terpart of the 7-ray one, with a delay of ~5 months due 
to opacity effects and the propagation of the shock along 
the jet. Past flaring episodes from PMN J0948+0022 have 
shown a different delay or no indication of correlated vari- 
ability llFoschini et al.l [201 ll : ID’Ammando et al .1120141 ). thus 
making the radio-gamma connection complex in this object. 
A slightly higher polarized emission from the core and po- 
larization percentage was observed by MOJAVE after the 
peak of the radio emission, while no significant change of 
the EVPA was detected. 

Following the 7-ray flare observed by the LAT, a 
first detection of VHE emission from PMN J0948+0022 
was attempted by VERITAS. The observations carried 
out between 2013 January 6 and 17, just at the end of 
the 7-ray flaring activity, resulted in an upper limit of 
-F>o. 2 TeV < 4.0 X 10 -12 phcm -2 s _1 . Further observations 
with the current generation of IACTs (VERITAS, MAGIC, 
and H.E.S.S.), especially during the peak of 7-ray activ- 
ity, may lead to VHE detections of this object. With a 
lower energy threshold and better sensitivity, future obser- 
vations with the Cherenkov Telescope Array will constrain 
the le vel of 7-ray emission at 100 GeV and below dSol et al] 
l2013h . However, pair production due to the EBL signifi- 
cantly suppresses the VHE flux of PMN J0948+0022 given 
its redshift of 0.5846. 1H 0323+342 is another, closer ( z = 
0.061), NLSyl, already observed (but not significantly de- 
tected) at VHE during a search for TeV-emitting FSRQs 
llFalcone et al.l [20041 ). If 1H 0323+342 exhibits flaring be- 
havior similar to that observe d in PMN J0948+0022 then 
it will be a promising target for IACTs. 

We compared the broad-band SED of the 2013 flar- 
ing activity state with that from an intermediate activity 
state of PMN J0948+0022 observed in 2011. Contrary to 
what was observed for so me FSRQs (e.g., PKS 0537—441; 
iD’Ammando et~al] l2013b|j the two SEDs, modelled as an 
EC component of seed photons from a dust torus, could not 
be modeled by changing only the electron distribution pa- 
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Table 5. Model parameters for the SEDs shown in Fig. [T] 


Parameter 

Symbol 

2011 intermediate 

2013 flare (dust torus) 

2013 flare (BLR) 

Redshift 

z 

0.585 

0.585 

0.585 

Bulk Lorentz Factor 

r 

30 

30 

12 

Doppler factor 


30 

30 

12 

Magnetic Field (G) 

B 

0.5 

1.0 

1.0 

Variability Time- Scale (s) 

tv 

10 s 

10 s 

10 s 

Comoving radius of blob (cm) 

K 

5.7xl0 16 

5.7xl0 16 

2.3 x 10 16 

Low-Energy Electron Spectral Index 

pi 

2.5 

2.35 

2.35 

High-Energy Electron Spectral Index 

P2 

4.0 

4.0 

3.9 

Minimum Electron Lorentz Factor 


2.7 

1.0 

1.0 

Break Electron Lorentz Factor 


5.0 x 10 2 

6.8 x 10 2 

7.8 x 10 2 

Maximum Electron Lorentz Factor 

'fmax 

1.0 x 10 4 

1.0 x 10 4 

1.0 x 10 4 

Black hole Mass (Mq) 

Mbh 

1.5 x 10 8 

1.5 x 10 8 

1.5 x 10 8 

Disc luminosity (erg s — 1 ) 

Ldisc 

5.7 x 10 45 

5.7 x 10 45 

5.7 x 10 45 

Inner disc radius (R g ) 

Rin 

6.0 

6.0 

6.0 

Accretion efficiency 

Vdisc 

1/12 

1/12 

1/12 

Fraction of photons scattered by corona 

fsc 

0.35 

0.35 

0.35 

Photon index of scattered coronal emission 

r sc 

2.3 

2.3 

2.3 

Seed photon source energy density (erg cm -3 ) 

'd'seed 

2.0 x 10 -4 

2.0 x 10“ 4 

0.3 

Seed photon source photon energy (m e c 2 ) 

€ seed 

7.0 x lO" 7 

7.0 x 10- 7 

3.7 x 10“ 6 

Dust Torus luminosity (erg s — 1 ) 

Ldust 

1.1 x 10 45 

1.1 x 10 45 

1.1 x 10 45 

Dust Torus radius (cm) 

Rdust 

3.8 x 10 18 

3.8 x 10 18 

3.8 x 10 18 

Dust temperature (K) 

- L dust 

1700 

1700 

1700 

Jet Power in Magnetic Field (erg s — 1 ) 

Pj.B 

5.4 x 10 4S 

2.2 x 10 46 

5.5 x 10 46 

Jet Power in Electrons (erg s — 1 ) 

Pj,e 

6.3 x 10 44 

2.4 x 10 45 

1.5 x 10 44 


rameters. A higher magnetic field is needed for modelling 
the high activity state of PMN J0948 +0022, as wa s also ob- 
served for the FSRQ PKS 2142—75 dDutka et all I2013l 4 . In 
a Compton dominance vs. synchrotron peak frequency plot 
(see Figure [5]), the values for the two activity states of PMN 
J0948+0022 lie in the same region occupied by the FSRQs, 
as well as the other NLSyl SBS 0846+513. This suggests 
that the EC scattering is likely the dominant mechanism for 
producing high-energy emission in PMN J0948+0022, con- 
firming the similarities between 7 -ray emitting NLSyls and 
FSRQs. The inverse Compton (IC) peak of this source seems 
to lie at 100 MeV or below. With the improve d sensitivity 
of the LAT at low energies with Pass 8 data (iAtwood et ali 
l2013h we will be able to characterize in more detail the IC 
bump of PMN J0948+0022, and therefore the SED mod- 
elling. 


In addition to the large uncertainties on the BH mass o f 
the NLSyl ('e.g.. lMarconi et ahll2008l : [Calderone et al.l[2013h . 
the structure of the host galaxies is unclear. For four out of 
five NLSyls detected in 7 rays by LAT, no high-resolution 
observations of their host galaxies are available. This is also 
related to the redshifts of these objects, and then the dif- 
ficulty of detecting significant res olve d structures for their 
host galaxies (e.g., SBS 0846+513: iMaoz et al.lll993j ). There- 
fore we cannot rule out the possibility that these 7 -ray emit- 
ting NLSyls are hosted in elliptical galaxies. 
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